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The Hamiltonian describing the system under discus- 



We develop a scheme to exactly evaluate the correlation 
energy in the random-phase approximation, based on linear 
response theory It is demonstrated that our formula is 
completely equivalent to a contour integral representation re- 
cently proposed in Ref. Q being numerically more efficient 
for realistic calculations. Numerical examples are presented 
for pairing correlations in rapidly rotating nuclei. 



Mean field theory provides a powerful approximation 
for describing many-particle systems. However, there are 
many situations where one is forced to go beyond mean 
field to include higher-order correlations, the random- 
phase approximation (RPA) being a widely used method 
for this purpose. One of the basic observables which re- 
quire to go beyond mean field approximation is the cor- 
relation energy. In this case, the different RPA modes 
contribute in a democratic way. To calculate this quan- 
tity one has, therefore, to determine many RPA eigen- 
modes, especially in the case where symmetries of the 
mean-field are spontaneously broken, such as in the case 
of superfluid and deformed rotating nuclei. In Ref. Q we 
developed a method to calculate the correlation energy 
in the RPA, making use of the response function tech- 
niques, and applied it to the study of pairing correlations 
in rapidly rotating nuclei. The essence of the method con- 
sists in expressing the RPA correlation as an integral in 
terms of the RPA response function, function which can 
be calculated without explicitly solving the RPA eigen- 
value problem. These techniques have been recently ex- 
tended to deal with the Nambu-Goldstone modes ||, and 
to calculate the nuclcon effective mass in superfluid, de- 
formed, and rotating nuclei || . 

Recently, a contour integral representation for the RPA 
correlation energy has been proposed in Ref. It was 
claimed that this method is more efficient than that de- 
veloped in Ref. Q. The assertion was also made that 
the method of Ref. Q did not take the contribution of 
the spurious modes into account. In this letter we show 
that both methods are completely equivalent. Further- 
more, using the property of meromorphic functions, we 
improve the original method of Ref. |j| to obtain the ex- 
act RPA correlation energy in more efficient way than 
that proposed in B . 



H = H + V, 



(1) 



where Hq is the unperturbed one-body (mean-field) 
Hamiltonian and V is the residual two-body interaction, 
which is assumed to be of the multi-separable form 



V = - ^^XpQpQp, 



(2) 



Q p being a one-body hermitian operator while Xp is the 
strength of the interaction in channel p. The associated 
ground state energies and state vectors of H and H are 
denoted E , |$ ) an d E, \^5>), respectively. Turning on 
the interaction adiabatically, the correlation energy can 
be written as || 



E 



E - En 
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(3) 



In this equation l'F(A)) is the ground state of the A-scaled 
Hamiltonian H(X) = Hq + XV. Within the RPA approxi- 
mation, the correlation energy takes the form (see e.g. |J) 
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where u> n is the RPA eigenfrequency and E a + Ep is the 
unperturbed two-quasiparticle energy (eigenstates of H a ) 
in the quasiparticle representation. As shown in H, the 
above expression can be rewritten as 



E, 



where 
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F(w)= / d\Ti-\ll w (uj)x 



(5) 



(G) 



The A-scaled RPA response function (matrix) TZ^'(u)) 
is defined in term of the unperturbed response function 
(matrix), 
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FIG. 1. An illustration of the integration contour in the 
complex plane. Crosses denote the positions of the RPA and 
unperturbed roots, which give singularities of the funtcion 
F(z), Eq. (|). 



as 



R (X) (w)= [l-iJHxA]- 1 ^), 



(8) 



where q p (af3) = (a/3\Q p \0) and \ = (<WX/=>)- 

A small but finite value of the imaginary part e = 
Im(cj) has been used to evaluate the correlation energy 
in Ref. (cf. Eq. ||) and found that it takes increasingly 
more computational time to approach the exact value of 
^corr^ ( e ^ 0)- Here we show that Eqs. @— (||) are equiv- 
alent to the integral representation proposed in Ref. m , 
and present a more efficient way to evalute the associ- 
ated correlation energy than that presented in Ref. 0]. 
For this purpose, we make use of the analytic properties 
of the function F(u>), which become evident by carrying 
out the A-integration in Eq. (o) analytically, 



F(z) = -log(dct[l-i?(z) X ]). 



(9) 



Thus, F(lu) has logarithmic singularities at the RPA 
eigenenergies and at the two-quasiparticle unperturbed 
energies on the real axis. Eq. (||) can now be expressed 
as contour integral (see Fig. |l|) 
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dzF(z), 



(10) 



where CQj, is a (closed) path which passes through the 
origin of the complex plane z, and which encloses all the 
positive RPA and unperturbed roots. Even if there is 
a spurious zero-energy mode (the symmetry-recovering 
mode or the Nambu-Goldstone mode) in the RPA spec- 
trum, F(oj) ~ const, x log z at z ~ so that the integral 
converges at the origin. In order to show the equivalence 
of Eq. ( JTo| ) to the corresponding formula of Ref. Q , we 
consider the integration path dipicted in Fig. [l]b, where 
a is a small positive quantity chosen to be smaller than 
the lowest singular point of F(z). Integrating by parts, 
one obtain 
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F(z). (11) 



If the zero mode exists, the segment of the real 
axis between the origin and the lowest RPA (or two- 
quasiparitcle) root is a "branch-cut" of the complex log- 
arithmic function F(z), and so [zf (z)] ^ = —2iria (mi- 
nus sign arizing from the fact that the direction of the 
path is clockwise). On the other hand, the singularities 
of the function z-4-F(z) are poles at the same points as 
those of F(z). Thus, in the limit a — > 0, the branch-cut 
contribution vanishes, and 



^ = ^{ dzz4-F(z), 



dz 



(12) 



which is nothing else than the integral representation of 
£^ r A of Ref. 0] (cf. Eq. (6) in it). Note that if there is 
a zero mode, ^F(z) ~ 1/z 2 , so the integral in Eq. ( |T2| ) 
diverges when the path goes through the origin. Other- 
wise, the origin is not a singular point and the path C|i|b 
can be trivially modified into in Eq. ( |l2| ) . If one uses 
the path Cjjjb in Eq. ( |To| ) instead of C^, one has to add 
the branch-cut contribution —a/2 if the zero mode exists. 

In spite of the equivalence mentioned above, Eq. (|Io| ) 
is numerically easier to calculate than Eq. (jl^) (or equiv- 
alently Eq. (6) of Ref. [Q), since the meromorphic func- 
tion F(z) has a more regular asymptotic behavior than 
z£F(z). In fact, in the limit \z\ -> 00, F(z) -> o(l/z 2 ) 
while z-j^F(z) — » o(l). Consequently the contribution 
to the corresponding integral arising from the segment 
A2B2 in Fig. vanishes when the points A2 and B2 are 
taken to infinity. Now the meaning of the approximation 
used in Ref. 0, is clear; the calculation of Eq. (||) with 
finite e = Im(cu) is equivalent to the integration along 
the path shown in Fig. ^a), except the contribution from 
the segment A1B1 which vanishes only in the limit of 
e = Im(w) — > 0. We have checked this point numerically. 




FIG. 2. The modified integration contour from Fig. 
which is more suitable to calculate the RPA correlation en- 
ergy. One has to take the limit of infinite radius of the semi- 
circle, whose contribution then vanishes. 

For convenience, the integration path in Eq. (|l^) can 
be modified to the one shown in Fig. ||a). In this case 
the contribution from the semicircle vanishes as its ra- 
dius goes to infinity. Moreover, the modification of the 
path parallel to the real axis into the one parallel to the 
imaginary axis is very useful for making the numerical 
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calculations efficient. This is because ImF(z) is an oscil- 
lating function of Re(z) on the path shown in Fig. [j]a), 
while ReP(z) is a rapidly decreasing function of Im(z) on 
the path shown in Fig. ||a). Consequently, one needs in 
this case a smaller set of mesh points to carry out the in- 
tegration than in the case of the path shown in Fig. |l|a). 
Therefore, the expression 
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duRe[F(iu))], 



(13) 



of the RPA correlation energy is numerically more con- 
venient to evaluate than the expression displayed in 
Eq. © or (|l|), equivalent to Eq. (6) of Ref. {§. It 
is worthwhile noticing that the integral appearing in 
Eq. ( |Io| ) or (|l3|), is further simplified by using the 
following properties of F(z): [F(z)\* = F(—z*) and 
F{—z) = F{z), which can be easily demonstrated mak- 
ing use of Eqs. (|J) and ([)[). Consequently the integrand 
ReP(z) is symmetric with respect to the real axis, and 
— 2 J °° in Eq. (jl3|). A similar simplification is pos- 
sible in Eq. (jlO|). 

In Refs. the pairing correlations in rapidly rotat- 
ing nuclei have been studied by using the general method 
discussed above. In these references, in addition to the 
RPA correlation energy, another measure of pairing cor- 
relations was introduced, namely the RPA pairing gap, 
^rpa 01 (called the "effective" pairing gap in Q). It is 
defined as 



-iia-s = \M 2 + ^ 2 S (RPA), 



with 



So(RPA)^ Yl [l(n|P|0)| 2 + |(n|pt| )p 



RPA 



(14) 



(15) 



where A = G(0|P t |0) H B is the standard, static BCS 
pairing gap (the order parameter of mean-field) , while G 
is the pairing force strength. The non-energy weighted 
sum rule So (RPA) describes the contribution of the RPA 
fluctuations for the monopole pair transfer operator, 
= J2i>o c l c l- Note that X)«^ng means that the di- 
vergent contribution from the spurious mode (pairing ro- 
tation) is to be excluded, in keeping with the fact that its 
contribution to Eq. |l4 is included through the static pair- 
ing gap A. In Ref. |1], So (RPA) was calculated making 
use of the expression 



S (RPA) 



dui Im Tr \TZ(uj + iej] , 



(16) 



where K(ui) = TZ^ =1 \uj) is the RPA response function, 
whose dimension is 2 corresponding to Qi = (Pt+P)/v2 
and Qi — i(P> — P)/>/2. A finite value of e and a low- 
energy cutoff w cu t are used to get rid of the NG mode con- 
tribution numerically. This is the same approximation as 



that used in calculating the RPA correlation energy 
and can then be avoided in keeping with the discussion 
leading to Eqs. © and (|l|). In this case, the path in 
Fig. ||b) is to be used in order to avoid the singularity 
associated with an eventual zero mode, as in this case 
TZ(z) has a second order pole at the origin (cf. ||): 

1 r + °° 

S (RPA) = — / dojReTr[U(a + iuj)]. (17) 

2 71 " J-oo 

Since the function Tr [72.(^)1 has poles as singularities, 
the integral is independent of the choice of a. Summing 
up, making use of Eqs. flU) and @, both the RPA cor- 
relation energy and the RPA pairing gap can be exactly 
evaluated in a numerically efficient way. 
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FIG. 3. RPA Pairing gap (upper panel) and RPA correla- 
tion energy (lower panel) for neutrons in 164 Er as a function 
of the rotational frequency. Both quantities are in MeV. The 
solid, dashed and dotted lines denote the results of calcula- 
tions with e = Im(u) = (exact), 100 keV and 200 keV, re- 
spectively, ftwcut = 400 keV is used for the approximate (finite 
e) calculations of the RPA paring gap. The value of the static 
(mean-field) pairing gap A, which vanishes at hu} Iot = 0.34 
MeV, is also displayed in the upper panel. 

In what follows we compare the results of the exact and 
approximate calculations of both an d ^RPA in the 

case of deformed, superfluid nuclei as a function of the ro- 
tational frequency (cf. Fig.(3)). Eqs. (HI) and (0) were 
used to carry out the exact calculation, while Eqs. (0) 
and d||) with finite values of e = Im(w) = 100, 200 keV 
were used to obtain the approximate results, as was been 
done in Refs. (Im(w) = 80 keV was used in |],0). 
The nucleus 164 Er has been chosen as a typical rotating 
nucleus and constant deformation parameters were used 
for simplicity. Only the monopole pairing force has been 
included with the smoothed pairing gap method being 
employed, which leads to a slightly different model space 
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from that used in JjQ]. Note that the correlation energy 
given in Eq. (^) as well as the sum rule value of the pair- 
ing gap given in Eq. ([l5]) include the exchange (Fock) 
contribution. In Fig. [| this contribution is excluded for 
the RPA pairing gap in accordance with Ref. J?j , while it 
is included for the RPA correlation energy as in Ref. 0] . 
The cusp behaviours are clearly visible in the correlation 
energy, which are caused by the sudden change of the 
mean- field (static pairing gap); they are associated with 
the pairing phase-transition at fi.w ro t = 0.34 MeV, and 
with the crossing of the g- and s-bands at 7iu; ro t = 0.25 
MeV. They are more evident in the RPA pairing gap, 
where the "exact" (e = 0) calculation diverges at these 
transition points. These divergences are due to a contri- 
bution of the lowest solution, which approaches to zero- 
energy at the transition points and brings about similar 
effects as those caused by the spurious mode, e.g. infinite 
strength. This is the well-known drawback of the RPA, 
whose small amplitude approximation breaks down near 
the transition points. The approximate results (e ^ 0) 
are very similar to the exact one (e = 0) and provide ac- 
curate estimates of the rotational frequency dependence 
of both the correlation energy and the pairing gap, ex- 
cept at the transition points. In particular in the case of 
the pairing gap, the singular behaviour at these points 
are smooth out because of the low-energy cutoff ?ko> C ut- 
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FIG. 4. Comparison of the RPA and the number projec- 
tion (NP) calculations of the pairing gap (upper panel) and 
correlation energy (lower panel) for neutrons in 164 Er as a 
function of the rotational frequency. The solid (dotted) lines 
denote the results of the RPA (NP) calculations. The results 
of approximate e = 200 keV calculation in Fig. [] are also 
included as dashed lines. 



compared in Ref. , where the RPA calculation was car- 
ried out approximately as in [Q. Here we compare the 
NP results with the exact RPA results in Fig. |. The 
NP correlation energy is defined as the energy difference 
between the NP and mean-field (Hartree-Bogoliubov) , 
Eco^ = -Enp — -E-HB (the exchange energy is included in 
E^p). Although RPA leads to larger values of the corre- 
lations, especially in the superfluid phase, the rotational 
frequency dependences are quite similar as already found 
in JjJ . The advantage of the NP method over the RPA is 
to lead to smooth functions for both the correlation en- 
ergy and the pairing gap at the pairing phase-transition 
point. It is, however, noticed that the cusp behaviour re- 
mains in the correlation energy at the g-s crossing point. 

In conclusion, a method to exactly deal with RPA cor- 
relations, based on that proposed in Ref. has been 
developed. It is equivalent to a recently proposed inte- 
gral representation in Q , and allows one to calculate the 
exact RPA correlation energy and in a numerically effi- 
cient way, properly dealing with the contribution of the 
Nambu-Goldstone modes. Making use of this method the 
pairing correlation in rapidly rotating nuclei can be stud- 
ied in detail, providing the basis for an eventual analysis 
of the pairing phase-transition in strongly rotating nuclei. 
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There is an another method which allows to go beyond 
mean-field approximation, namely number-projection 
(NP) (see e.g. ||). The results of both methods were 
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